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MATHEMATICAL BIOPHYSICS OF GROWTH 


N. RASHEVSKY 
The University of Chicago. 


The rate of growth of a tissue is studied mathematically in its 
dependence on the metabolism of the cells. A high glycolytic coeffi- 
cient, which facilitates cell division, as has been shown before, does in 
this way also increase indirectly the rate of growth of the tissue. 
There is however also a possible direct effect of glycolysis on the 
rate of growth, which is also studied analytically. Equations are de- 
rived, giving the total rate of growth of a tissue in its dependence on 
the glycolytic coefficient. 


In our book (Rashevsky, 1938, hereinafter referred to as MB) we 
have derived some relations (MB, Chapt. X), showing that the forces 
which tend to divide a cell are the stronger the greater the glycolytic 
coefficient £ of the cell, and pointed out, that this is in a certain agree- 
ment with Warburg’s findings, that abnormally rapidly multiplying 
cancer cells exhibit an abnormally high glycolytic coefficient. An ab- 
normally easy cell division is however not directly connected with an 
abnormal rate of growth of the cellular aggregate as a whole, and it 
is the abnormal rate of growth that is especially characteristic of 
tumors. In the present paper we shall give some general theoretical 
discussions of possible relations between glycolysis and rate of growth. 


P L 

4 Following the method developed in the Appendix of MB., we shall 
consider cells of an arbitrary shape. When they are not too elongated, 
- they are characterized by their average overall size a. 

Let a substance enter such a cell and be consumed there at a rate 


Q proportional to the concentration c, 
: Q=ke. (1) 
Denoting by ¢ the external concentration, by c the average internal 


concentration, by ¢, the internal concentration at the membrane, and 
using the same argument as before, we have (Rashevsky, 1939) : 


Anha? (Co — C1) = fs natke (2) 
and 
—119— , 
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—c 
daha (to — 0.) = Se —* ; (3) 
or 
ké 
h (CG) (4) 
3 
and - 
h( Se aa (5) 
Combining (4) and (5) and rearranging: 
ak .— 
= (1+ ge: (6) 
Substituting (6) into (4) we find: 
oa Co 
aH : wk ak (7) 
6D! 3h 


If the molecules of the substance are rather large (amino-acids), 
so that D and h have small values, of the order of D ~ 10-*° cm? - sec™; 
he 10-*cm.sec?, we find fora~ 10° cm, k ~ 10°sec* ; ¢ ~ 10% 
gr-cm= 


azk ak 
Ait Chee oe (8) 


For simplicity consider the special case, in which D ~ 10-1°cm?sec* 
but h >> 10-*cm - sec. Then equ. (7) reduces to: 


AW ates (9) 


Let the substance be transformed into another one, forming the 
body of the cell (proteins). Let the latter break down at a constant 
rate q gr.cm-*.sec. Denoting by V the volume of the cell, we then 
have for a constant concentration 6 of that other substance: 


dV = 


Introducing (9) into (10) and remembering that V = < att, we 
find: 


rere (11) 
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The size a tends to a limiting value a, obtained by making the right 
hand side of (11) equal to zero. This gives: 


6De, 
— 


ES (12) 


Taking for g~ 10° gr .cm®. sec", a value found for the breakdown 
of cell proteins (MB, p. 104) we find: 


aw 102em. (13) 


In actual cases not one, but several substances (amino-acids) are 
consumed, building up cell proteins. Let there be m substances con- 


sumed each at a rate k;c,c,---c,. For the i-th substance the rate of 
consumption may be written: 


ke; (G15 --- Cia Cisne? Cn) C3 ° (14) 
The expression i; (C12 --+ Gi Cis --- Cn) plays the role of k in equ. (1). 
By a similar argument as before we find approximately 
= 6D;c A 
Ce aon a (15) 
O71, Co +++ Cin Cis >>> Cn 
where ¢; is the external concentration of the i-th substance. Hence: 


are a 6"D,D, = D€01Co2 *** Con 


DE CL SSS SS (16) 
A?" Hehe +++ Ken (C1Co +++ Cn) 
or 
ae 67/D,--- Dn VCor*** Con __ a (17) 
OF 4/ Ke, +++ Kn ak 


which is of the same form as (9). For brevity we shall write here 
_ D, k etc. as if one substance only were present. 

From the discussion in Chapt. X of MB, which leads to the 
inequality (8) on p. 105, it follows that the resulting “equivalent” 
outflow is of the form 


Ag ae (18) 


where f is the glycolytic coefficient and A and B are functions of 
temperature, diffusion coefficients and molecular weights. 
The critical stability size a* is given (MB, Appendix) by 


at act 9MD,y =n 9MD,y 
RTqu VN RT(AB—B)p 


(19) 
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where M and D, are some sort of averages of the values for the 
different metabolites involved in the respiratory metabolism. The 
values of D, is different from that of D, the latter referring to amino- 
acids. 


If a* >a , as given by (12) then the cell does not divide at all. 


Ti Ga a, then a cell divides as soon as a = a’. 
Equ. (11) gives: 

2 

2De— 5 a? 


which after elementary integration and rearrangement gives: 


2q 
-—t 
OD Com. i, 36 , (20) 


a? 
where C is an integration constant. 
If the cell divides at a = a*, then each half cell has the size 
a = 0,795 a*, or 
G7 0,63 625 (21) 


At this size it starts growing again. We take the moment at which 
a=a ast =O. Hence fort = 0 


D 
Ene Wess. 
q 
or 
6D 
C= = — 0,63 a". (22) 
Hence, introducing (22) into (20): 
2q 
6De 6De ae 
Bis ee (ee Soca eae (23) 


The value At of the time interval between two successive divisions is 
found by equating the right hand side of (23) to a® (since for t = 0 : 
a =a’), and putting t = At. This gives after some rearrangements: 


2q 
-—At 


tine oe 6De, — qa*? 


6Dc, — 0,63qa*2 ’ 
or 
36 y 6Dc, — 0,63 ga” 
2q £ 6De,2=7a8e U (24) 
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If in an aggregate of N cells every cell divides every At seconds, then 
from N cells every At seconds N new cells are made. Hence 


dN WN 
iby winds 
or 
t 
N=N,e. (25) 


Introducing (19) into (24) and the latter into (25), we obtain 
the rate of growth of a cellular aggregate as a function of the gly- 
colytic coefficient 8. The larger £, the smaller a* , the smaller Aé and 
the more rapid the growth. The increase of a with time between divi- 
sions is practically linear for large 6’s; convex upwards for small f’s. 

Equ. (24) may be written: 


0.63 qa‘? 
36 NII 
AE Ses Age Te See: 
reeryare 


With values used above, ga*?/6Dce,.~ 10° << 1. Hence (26) may be 
written: 

_ 0.09 6 
Deg 


Taking 6~ 10", as representing the order of magnitude of the dry 
weight of a cell, we find with the values used here for the different 
constants, At ~ 1 — 10 days. 

Introducing (19) into (27) and the latter into (25) we find: 
Deo 4 RT(AB-B) \22 
N= Neewws (aay) : (28) 
In this simple theory, a tissue grows the faster, the greater /, other 
conditions being equal. But a faster growing tissue should have also 
smaller cells, other conditions being equal. It may be remarked that 
cells of an early embryo are usually smaller than cells of an adult or- 
ganism. 


At a*? . (27) 


II 
Glycolysis may however influence the rate of growth in a more 
direct way, than that considered above. When f is large, the net ef- 
fective metabolic flow is directed outwards and results in a pressure 
on the cell membrane directed outward. For an approximately con- 


stant q, this pressure is equal (MB, Chapt. VIIZ) to 
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(29) 


where h again refers to an average value for respiratory metabolites. 
The pressure (29) may be partially compensated by a static osmotic 
pressure ,, due to different nonmetabolised substances, dissolved in 
the cell and in the surrounding, for which the membrane may be im- 

permeable. The total pressure then is equal to 
bh lake (30) 

M 3h 
Let a reversible reaction take place in the cell: 

A, +A,+---AnaP,+ P.+---+ Pa, (31) 


where A; represents some aminoacid, and P; a protein. If this reac- 
tion requires a special enzyme, which is present only inside of the cell, 
then it will not take place outside of the cell. If the membrane is per- 
meable to A; but not to P; , then the concentrations of A; will be the 
same inside of the cell and outside, and if the cell size remains con- 
stant, the amounts of P; will not change, everything remaining static. 
Let however, due to some cause, the cell size increase very slowly. 
Then the concentrations of A; and P; will decrease. But due to the 
supply from outside, the concentration of A; will be practically imme- 
diately reestablished. Thus the concentration of A; will remain the 
same, that of P; decrease. As a result of this some P; will be formed 
according to (31). The cause of the increase may be the pressure 
(30), which tends to expand the cell. According to this picture a cell 
does not grow in size because of increase in mass, but rather increases 
in mass because of growth in size. 

In the first approximation the increase of size a will be propor- 
tional to (30) and given by 


=~ =o (=o) (32) 


where ¢ is principally determined by the resistance of the membrane 
to the flow of water, which due to its incompressibility must flow into 
the cell as soon as the latter expands. Equ. (32) integrated gives: 


Br ee se 


C being an integration constant. 
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For 6 just above the critical value 6. = B/A , when q = AB —B 
is positive but very small, such a cell might grow to tremendous di- 
mensions, though the rate of growth would be extremely slow. For 
normal f it will grow to a normal size a and then divide. Again, as 
before, let us put ¢ = 0, when a = 0, 8a*. Hence 


3Mhp, 


= Ci 058.0", 4 
RTq mi a) 
or 
Gowan oe (35) 
RTq 
Therefore 
ss 2 RT eq 
3Mhp, Mhp, —* 
c= —+ (0,8 a* — ues? ae (36) 
Tq q 
The interval At between two successive divisions is given by 
a= Cow RT eq 
38Mhp, SV Rese ao 
—-| (0,8 a*— ——_—__) e=™* =a’, 
RTq RT 
or solving for At 
3Mh RTqa* — 3Mh 
= Sea Cie See oe 2 . (37) 
RTlq 0.8 RTqa* — 3Mhp, 
For very small values of p, , (37) reduces to 
oe Vi 
Ate ey log 1,25 es iGTEET ; (38) 
RT¢q BICq 
which introduced into (25), gives: 
RT¢q : 
N= N, e 0,67uh . (39) 


As p» increases and approaches the value 0.38RTqa*/3Mh, At ap- 
proaches infinity, because for large values of p. the net pressure (30) 
is directed inward and the cell does not grow. 

The larger f, the larger q and the larger the rate of growth. The 
cell size a* does not need however to be smaller in this case, in spite of 
a larger q, since (39) does not contain D,. If D, increases in the 
same ratio as g/M , a* remains unchanged. This holds a fortiori if in- 
stead of (38) we use the more exact expression (37). 
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An increase of ¢ in (39) results in an increased rate of growth. 
On general physical grounds we would expect h and C to increase si- 
multaneously, though not necessarily in the same ratio. Nothing can 
be said therefore about the influence of h on the rate of growth, until 
we possess a better developed theory of permeability. From what was 
said above about the role of D, , it follows that more rapidly growing 
cells should have a larger D, , if their size is not smaller than that of 
slower growing cells. Inasmuch as cancer cells are found of any size 
within normal range, we may expect them to possess a higher D, 
which would agree with the observations concerning their higher 
water content and higher general permeability. 


III 
We may also consider other limiting factors for the cell size than 
the mechanical instability or the one given by equ. (12). Let the res- 
piratory reaction be catalized by some enzyme, whose concentration is 
m gr-cm=, so that 


q=am, (40) 


a being a coefficient of proportionality. Let m vary according to the 
equation: 


ee Klett g's (41) 


where q’ is a constant and c’ is the concentration of a substance, which 
flows into the cell, and is there transformed at a rate k’c forming the 
catalist. In this case ¢ is given by (9). Hence 

dm _ 6D'c, 


“dt a? 
D’ and ¢,’ in (42) are of course different from D and ¢ in (9). Equ. 
(42) will hold with good approximation even when the size of the cell 
is not constant, but increases with time very slowly, so that the quasi- 


stationary diffusion equation may be applied (MB, Chapt. I). 
For 


—q. (42) 


6D'c,’ 


rien 


or 


6D'c,’ 
a>,/ Fs Me (43) 
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dm/dt becomes negative. Hence m and therefore also gq = am will 
decrease, when the cell exceeds the size m. For sufficiently large a , q 
will rapidly become so small, that RTqa/3Mh = p,, and therefore 


da/dt = 0, because of (32). 


Or m may vary like 
dm = : 
—— = Ko’ —k"m. (44) 


Again if a does not vary too rapidly, we have for each a the quasi-sta- 


tionary value: 
6D'c,’ 
a (45) 


= = 
ka? 


Hence from (45) and (40) 
6aD'c,’ 
‘ (46) 


Because of (32), the cell will not grow above the size, given by 


RTqa 


3Mh 


ae BOs (47) 


Introducing (46) into (47) we find for the limiting size 
2RTaD'c,’ 
se Marieke (48) 
The author is indebted to Mr. H. D. Landahl for the checking of 


calculations and for several critical remarks. 
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STUDIES IN THE MATHEMATICAL THEORY OF EXCITATION 


A. S. HOUSEHOLDER 
The University of Chicago 


The general linear two-factor nerve-excitation theory of the type 
of Rashevsky and Hill is discussed and normal forms are derived. 
It is shown that in some cases these equations are not reducible to 
the Rashevsky form. Most notable is the case in which the solutions 
are damped periodic functions. It is shown that in this case one or 
more—in some cases infinitely many—discharges are predictable, 
following the application of a constant stimulus S. The number’ of 
discharges increases with S , but the frequency is a constant, charac- 
teristic of the fiber and independent of S. 


1. The general linear two-factor theory. The two-factor nerve- 
excitation theories of Rashevsky (1933) and of Hill (1936) are nat- 
ural generalizations of the single-factor theory of Blair (1932), which 
supply several of the deficiencies of Blair’s theory. At the same time 
these have not appeared to yield repetitive discharges of the nerve 
fiber under constant stimulation, a phenomenon which is often met 
with empirically. It is the purpose of the present discussion to con- 
sider the most general possible linear two-factor theory, and to show 
in particular that for suitable choices of the parameters of the equa- 
tions such repetitive discharges are predictable. 

In its most general terms the two-factor theory postulates the ca- 
pacity of the nerve fiber to develop two “substances” or “factors”, the 
rate of development of each being a linear homogeneous function of 
the three quantities: the stimulus intensity, and the excess of each 
substance or factor over the resting value. Excitation is supposed to 
occur and to continue as long as a certain linear homogeneous function 
of the measures of these factors is positive. 

We shall speak of the factors as substances, for convenience and 
definiteness of the picture, though they do not need to be such, and 
we shall speak of their concentrations as measures of the factors. 
Then if z, and x, are the concentrations at any time ¢, and if S(t) 
is the stimulus intensity, then the development of the substances is 
governed by the linear differential equations 


—129— 
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Ee = yi (41 a X1°) s\ Ay2 (X2 — 2°) mf Oy S(t), 
(1) 
sa = Oo, (%, — 04°) + Ao (Lz — H2°) + Oe S(t), 


where x,° and z,° are the concentrations in the resting fiber. Since we 
are by no means insisting that they really are substances being devel- 
oped, we shall not require that x, and «2 be positive to be meaningful. 

By suitable choice of units and of subscripts, it is no restriction 
to assume that the condition for excitation be of the form 


L,—%>O0, (2) 


the left member of the inequality being the linear homogeneous func- 
tion referred to above. 

Empirically only S(t) and the resulting interval of excitation are 
measurable, that is to say, only S(t) and the times at which the in- 
equality (2) is satisfied. Hence we shall define two forms of the two- 
factor theory as being equivalent in case the corresponding inequal- 
ities (2) are simultaneously satisfied. With this definition of equiva- 
lence we shall investigate the conditions for equivalence of any two 
two-factor theories and deduce normal forms for these. 

Blair’s theory is obtainable by setting di. = Ae: = Asoo = % = 0 
in (1). Rashevsky’s theory assumed adi2 = a2, = 0, while Hill’s the- 
ory had a1. = a. = 0. Offner (1937), seeking to test Rashevsky’s and 
Hill’s theory experimentally, found that they were, in fact, equivalent 
in the sense defined above. Young (1937) then showed that the most 
general two-factor theory (1) could in general be formally reduced to 
the Rashevsky form. However, when the characteristic roots are com- 
plex so are the resulting coefficients in the Rashevsky form. Physically 
this is the case of (damped) periodicity with possible repetitive dis- 
charge, and is most conveniently studied by reducing to a non-Rash- 
evsky form with real coefficients. Rashevsky (1938) has reviewed 
completely the case of the equations in his form with real characteris- 
tic roots, summarizing the equivalence proofs of Offner and of Young, 
and discussing empirical checks. 

2. Roots real and distinct. Equations (1) can be written in ma- 
trix notation in the form 


a 
FE U(a— 2°) +a S(t), (3) 


where S is a scalar, a, « and x are column vectors, and a is a two-by- 
two matrix. Any linear substitution 
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y=cu, —=c—, (4) 


where the matrix c is a non-singular matrix of constants, transforms 
the linear differential equations (3) into the linear differential equa- 
tions 


d 

= b(y—v") +B S(t) (5) 
where 

pe="CC*s b=ca. (6) 


However, we can admit only those matrices c for which inequality (2) 
and the inequality 


Y¥i1—Y2>0 (7) 


are simultaneously satisfied. Such matrices c will be said to define an 
admissible substitution. Any equations (3) and (5) obtainable one 
from the other by an admissible substitution are equivalent in our 
sense. 
It is at once evident that a scalar matrix 


e=(5 ») See (8) 
p 

defines an admissible substitution. This has the effect only of multi- 

plying the two coefficients a, and a, by the same positive scalar factor 

p- Hence, only the ratio a, : is important, and we may, for example, 

at any time assume o. and a, to be the sine and the cosine of some 

angle. 

We next recall the well known theorem in algebra which states 
that for any non-singular matrix c, the characteristic roots of the ma- 
trices @ and c a c are the same. These are the roots 4, and A, of 
the quadratic equation 


tig diese (9) 


where I is the identity matrix. Hence if an admissible substitution 
exists such that the matrix b is diagonal, b has necessarily the form 


Ae. <0 
f a (10) 
We can easily write down, the matrix c defining such an admis- 
sible substitution, when the roots 4, and A, are real and distinct, by re- 


ferring to some principles of projective geometry. Consider the trans- 
formation 
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F&=aég (11) 


of the elements (é:, 2) of a one-dimensional projective form into the 
elements (,', 6’) of this same form. When the roots 4, and A, are 
real and distinct there are two real and distinct fixed elements, 1.e. 
two elements of the one-dimensional form, (é:', §") and (€,2,, €27) 
which are transformed into themselves by the transformation (11). 
These are given by the two pairs of dependent homogeneous equations 


hi ise aot s G= bs 2), (12) 


If we introduce new coordinates into this projective form by the 
coordinate substitution 


n=cé, n=cé, (13) 
then the transformation (11) is equivalent to 
n =COC*n. (14) 


Now a projective coordinate system in a one-dimensional form is 
fixed when the projective coordinates of three elements of the form 
are assigned. Let us, therefore, assign to the point & the y7-coordi- 
nates (1,0), to & the y-coordinates (0,1) and to (1,1) the n-coordi- 
nates (1,1). Evidently, then, the points (1,0) and (0,1) are the fixed 
points of the transformation (14), and therefore this takes the form 


ni =Aini. (15) 


Since the required coordinate substitution changes the coordi- 
nates of €1 and & into (1,0) and (0,1) respectively, and leaves the co- 
ordinates of (1,1) unchanged, it is easy to write down this substitu- 
tion explicitly in terms of the é,;/ by expressing the fact that the anhar- 
monic ratio of an arbitrary & with &', & and (1,1) is equal to the an- 
harmonic ratio of the corresponding 7 with (1,0), (0,1) and (1,1). If 
we write, then, x and y in place of & and y for the variable point we 
obtain the desired form of the substitution (4): 


iA &,? SC &,° Vo _ ee vy — &43 Xe 
“dh waneereh cae tidmmoia tts af 4 aie See a 


where p is an arbitrary constant. 

One thing remains to be determined. We observe that interchang- 
ing the notations & and & of the two fixed points, or, what comes to 
the same thing, interchanging the subscripts on the two characteristic 
roots 2, and 4, , has the effect of interchanging y, and y.. Whichever 
root is called 4, and whatever the sign of p, the equations x, = x, and 
Yi = Y2 will be simultaneously satisfied, but unless these are properly 
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associated the orders of the inequalities will be reversed. For definite- 
ness we require that 


p>O. (17) 


It is no restriction if we so choose the homogeneous coordinates of 
&1 and & that 


&'—é,'=1. (18) 
Then we have 
Yo Uy — pt (&,4 <= § 27) U4 = (é1 a £,?) 2} . 
But by subtracting equations (18) one from the other we find that 
&' — &.2 = §,7 — és, 
whence 
Y1 — Yo= p (E22 — Fo") (2; — 22) . (19) 
Thus in order that the inequalities (2) and (7) shall be simultaneous- 


ly satisfied the designations A, and J. must be assigned to the charac- 
teristic roots in such a way that 


&2— &,1 >0 (20) 


when the scalar factor p is chosen positive and the homogeneous co- 
ordinates of the fixed points are chosen to satify (18). 
From equation (6) we have 
By = p(é.? Oh > he Oe) ’ 
(21) 
Bo = p(&2! Otero Eas Qs) . 
Varying the scalar factor p does not affect the matrix b , but only the 
magnitudes of the coefficients 6. Hence, we may choose this scalar 
factor so that 
Bi? + Bo? = (22) 
and hence so that f, and f, are the cosine and the sine of some angle. 
This is in accordance with the statement made above that only the 
ratio of the coefficients of S(t) is important. Note that 
Bb. — Be = p(&2? — §2*) (a, — Oe) ’ (23) 
so that the quantities 6, — fb. and a, — a have the same sign. Nothing 
can be said, however, about the signs of 6, and f, separately. . 
We have tacitly assumed, in the foregoing discussion, that the 
point (1,1) is not itself a double-point of the transformation (11). 
Postponing, for the moment, our consideration of this possibility, we 
may summarize: 
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Let the excitation equations in the explicit form (1) or the matrra 
form (3) be such that the roots A, and dz of the quadratic equation (9) 
are real and distinct. Let the vectors & and & , which satisfy the ma- 
trix equations (12), both have unequal components. Then their com- 
ponents can be chosen to satisfy (18) and the indices can be so ad- 
justed that (20) is satisfied. Then the substitution (16) with an ar- 
bitrary positive scalar p is adnussible and transforms the equations 
(1) into the equations 


eo =A(Y¥1 — ¥°) + Bi S (2) » 

(24’) 
Ys “fe: am 0)/. ’ 
a Ye )+ Bo S(t) 


where B, and B. are given by (21). The scalar p can further be speci- 
fied so that (22) is satisfied, and in this case the equations (24’') can 
be written 


se = 4,(Y1 — 1") —- S(é) cos B ? 
(24) 
= = dz (Ye — Y2°) += S(t) sin B ? 


for some angle 6. This is the Rashevsky form of the excitation equa- 
tions, and it contains five essential parameters. For stability to exist, 
and non-excitation in the resting state, it is necessary that 


Ay << 0 > dz << 0 ’ y,° Ss Y2° ’ (25) 


while from the nature of the substitution it follows that B, — B. and 
a, — ao satisfy (23) and hence have the same sign. 

In the exceptional case when (1,1) is a fixed point of the trans- 
formation (11) there is no substitution (4) admissible in our sense 
which throws the excitation equations into the Rashevsky form. If A, 
is the root corresponding to the fixed point (1,1) and if J, is the other 
root, the equations (1) are in the form 


i = @(x, — 44°) + (4, —@) (42 — 22°) + a, S(t) ; 


(26) 
a = (@— dz) (#1 — 1°) + (Ar + de — @) (22 — @,°) + 2 S(t) » 


where a is some constant. The other fixed point is then (a—dh, 
a — i4,) . We may make a substitution which gives to this point the 
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coordinates (0,1), in which case the excitation equations take the form 


dy, : 
aes Ai (Yi — yx) + B, S(t); 
(27) 
dY2 
dt = (A, Az) (Y1 — Yr") + A, (Yo — Y2°) + Bp, S(e) , 


where the coefficients 6 are yet to be determined and we may still 
assign the projective coordinates of another point. But £, and f are 
the projective y-coordinates of the point whose coordinates are 
(a, , a) in the é-system. Hence the following statement is immedi- 
ately evident: 

In the exceptional case when (1,1) is a fixed point of the trans- 
formation (11), the excitation equations cannot be given the Rash- 
evsky form by any admissible substitution, but they can be given the 
form (27). In these equations j, is the root of (9) corresponding to 
the fixed point (1,1), and A, is the other root. To obtain this form one 
has only to choose a substitution keeping fixed the coordinates of 
(1,1) and giving to the other fixed point the coordinates (0,1). As 
for the coefficients 6 , there are three possibilities: 

a. tf ao, = a, then BB, = B. whatever substitution of this type 
one employs; 

b. if (a,, a) is the second fixed point of the transformation 
(11) then B, = 0 and f. is — 1 or + lLaccording as a, — o& ts positive 
or negative; 

c. if the coefficients a do not satisfy either relation, then 6, and 
B. can be given arbitrary distinct values, 6B, ~ Bp. , and in particular 
fb. can be made equal to zero and f, equal to +- 1 or — 1 according as 
a, — a, is positive or negative. 

In case c the substitution is uniquely determined since the coor- 
dinates of the three points whose initial coordinates are (1,1) (a—d,, 
A>) and (a,, %2) are assigned. In cases a and b the point (a, , ae) 
coincides projectively with one of the other two, and the third point 
may be chosen at will. 

3. Roots real and equal. If the characteristic roots are equal, 
4, = 4s = A, and the matrix b can be diagonalized, then it is a scalar 
matrix 


Be EAL 
eact=b=( A) 


whence 
a=c1bc=b 
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since a scalar matrix is commutative with any matrix. Hence in the 
case of equal roots the excitation equations cannot be given the Rash- 
evsky form unless they are initially in this form, and any (oumies 
substitution leaves them in this form. If, in addition, 0, = a,, then 
the two quantities 7, — x,° and x, — %,° satisfy the same differential 
equation, a case which is obviously of no importance, In the light 
of the discussion of the preceding section it is therefore evident that: 

If the characterictic roots are equal, then the excitation equations 
cannot take the Rashevsky form unless they are initially of this form. 
In this case, however, it is physically necessary that o, a2, and the 
equations can be given the form 


dy. , 
ai =1(Yy— yr’) S(t), 
(28) 
dYs ; 
ai = 1(Y2— Y2") » 


where the sign before S(t) in the first equation is that of the quantity 
01 — Oe. 

In the contrary case, when the excitation equations have equal 
characteristic roots but are not in the Rashevsky form, there is only 
one fixed element of the transformation (11). This can be given, say, 
the projective coordinates (0,1) by an admissible substitution, pro- 
vided it is not the element (1,1). The original equations in this case 
can be written 


da, 0 
Ge (AY) (1 — 21°) — w(t 


22°) +a, S(t), 


(29) 
OES eGR 
Ge bm) + Aa) (@2 — m*) + a2 S(t) 5 


where J is the characteristic root and (u, 1) is the fixed element. If 
a, = a, the substitution 


ba v1 — Ue 2 X11 — 4 
Ce epee Ye pee eee (30) 


On —— O4 


where 


p (a2 & —a,) (1— wu) (a, —a) > 0, (31) 
yields the form 
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dy - 
Tm A(Y¥i— vv’) = S(t), 
(32) 
dy : 
SF = b(y:— yi") tA(Y2 — ¥2°) » 
where 
( = iy) (a —— Os a 
2 OL (33) 
Oe Oy 
In case a, = a, we obtain the form 
dy, _ ‘ 
at =A(yi— yr) + S(E) ; 
(34) 
dy. — : ’ 
RF, ee b(y1— Y,°) + 1(Y2— Yo )+S(t) , 


with a different expression for b. If u = 1 in (29), the fixed element 
is (1,1) and the substitution yielding the form (32) or (34) is no 
longer admissible. 

If the characteristic roots 4, = i. = A are equal, and the excita- 
tion equations are not in the Rashevsky form, they can be written in 
the form (29), where (u,1) is the single fixed element. If u~ 1 and 
a, ~ a, a substitution of the form (80) yields the equations (32). If 
a, = a, it is possible to obtain the form (34). But if uw = 1, neither 
of these forms is obtainable by an admissible substitution and one can 
only alter the form (29) by making the coefficients a, and a, equal to 
zero or unity. 

4. The characteristic roots complex. In this case the roots are 
necessarily distinct, being conjugate complex. Let these bed + ip. 
Then we choose the substitution which gives the coordinates (1 , —%) 
to the fixed point corresponding to the root 4 -++ i uw, and the coordi- 
nates (1,7) to the fixed point corresponding to the root 1 —iwu. The 
excitation equations then take the form 


oe =1(Y. — 9°) — w(Y2 — Yo) + Bi S(f) » 

(35) 
dy, 
i u(Y1 — 2°) + A (Yo— Ye?) Se aS CL) ; 


The explicit form of the equations of substitution can be derived as 
explained previously, but they are slightly more complicated and we 
do not write them here. The substitution is uniquely determined up 
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to a scale factor p. It is convenient for the further discussion to ad- 
just the time units so that 


VP+tw=1. (36) 
This can always be done by substituting 
7 = ty Pee (37) 


and writing the corresponding equations in r instead of ¢. Suppose 
this has been done, and let us then rename the variables, calling the 
time variable in the new units ¢ instead of + and the new coefficients 
and functions againd, u, 6, and S(t). Equation (36) is then satis- 
fied and by a further choice of the scale factor p we can suppose that 
pr be = 1. (38) 
Hence we may set 

A=cosy, w=siny, 
(39) 

fb, = cos Bf, fb. = sin B, 


and write our equations (35) in the form 


dy, 
dt 


= (¥:— yr") cos y — (Y2 — y2®) sin y + S(t) cos #, 


(40) 
ae = (y:—y,°) sin y + (Y2— Ye") cos y+ S(t) sin B- 


The remainder of our discussion will deal with the properties of 
the solutions of the equations (40), or of equations (35) with (36) 
and (38) holding. 
If S(t) = 0, the homogeneous equations have the solutions 
Y¥1— yr = e' (A cosut—B sin wt), 
(41) 


y2 = e'(Beoosut+Asinut), 


as one can readily verify, with A and B constant. In fact, A and B are 
the values of yi; — yi°, and ys — y.°, respectively, at the time t = 0. 
In the general case the solution still has this form, except that the 
A and B are no longer constant, but are functions of ¢ defined by 


A=A,+ [e~s(a) cos (ua — 8) dz, 


Y2 


(42) 
B=B,— [e*S(x) sin (ua — B) dev. 
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In these equations A, and B, are constant, and are the initial values 
of y; — y,° and y, — y2° respectively. 

We are especially interested in the case when S is a constant and 
A, = B, = 0, i.e. when a constant stimulus is applied to a resting 
fiber. In this case the quadratures (42) can be effected, and one ob- 
tains 


A = S[cos(S — y) —e*' cos(ut— B+ y)], 
: (43) 
B= S{[sin(f— y) + e*' sin(ut—f+y)]. 


Hence, on substituting these values of A and B into (41) we obtain 


Yi—Yyv = S[e cos (ut + B— y) —cos(6 — y)], 
(44) 
Y2— Yo = Sle’ sin(ut + 6— y) —sin (6 —y)], 
or, somewhat more explicitly, 
Y: — yi° = S[e*°*7 cos(¢ sin y + 6 — y) —cos(6— y)], 
45) 
Yo — Y2® = S[e*s¥ sin(t sin y + 6 — vy) —sin(6— y)]. 
Excitation will occur at the first, third, --- , roots of 
: Yi = Y2 
i.e. of 
e**[cos (ut + 6B — y) —sin(ut+ f—y)] 
(46) 


= cos (6 — y) — sin(f — y) ae 


and will continue until the second, fourth, --- , roots respectively. This 
equation can be simplified slightly, by means of a trigonometric iden- 
tity, to the form 


& cos (ut + B—y +7) = cos(B—y+7) ere , (47) 
and if we set 
p=B+3 (48) 
_ we have 
Yo — y.° 
pct Ba 0 mes 3 Meh yer gee 49 
et cos (ut -+ B’ — y) = cos(p’— y) + oe (49) 


For abbreviation set 
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f(t) =e! cos(ut+ p—y), 
(50) 


paces § 0 
$(S) =# Yi Swi 


rege 
and write (49) in the form 
f(t) =f(0) + ¢4(S) . (51) 


The following statement is now immediate: 

When the excitation equations (1) have complex roots i a iu 
there is a unique admissible substitution (4) reducing these equations 
to the form (35), with B, and 6, satisfying (38). By a further choice 
of time unit it is possible to obtain the somewhat simpler form (40). 
In these equations 4 = cosy <0, u=siny 40, and y2° > yi, but the 
parameters are otherwise independent and unrestricted as to sign. Hx- 
citation is assumed to occur and persist while y1 > Ye. 

The general solution of these equations is given by (41) and (42), 
but in the special case that a constant stimulus S is applied to a rest- 
ing fiber the quadratures can be effected and the solutions are given 
by (44). The intensity-time relations are then given by an equation 
of the form (51) where f and ¢ are defined by (50), excitation lasting 
while £(t) —f£(0) > 6(S). Since 6(S) > 0 for all S , while vanish- 
ing asymptotically in S , and since f(t) fluctuates periodically in sign 
and vanishes asymptotically in t , we have the following possibilities: 

If £(0) > 0, f(t) can exceed £(0) at most a finite number of 
times, and if £(0) ts an absolute maximum of f(t) for allt > 0, no 
excitation is possible for any S. If, however, f(t) exceeds £(0) for 
some t > 0, then for sufficiently large S at least one discharge will be 
possible, and the rheobase is given by the value of S for which 
o(S) = f(t.) — f(0), t, being the value of t at the first maximum 
of f(t). 

If £(0) = 0 the rheobase §, is given by (S:) = f(t), and if ta 
is the time of the n-th maximum, n discharges are obtainable by mak- 
ing S > Si where $(S,) = f(ti). The solution S, exists and is 
unique foranyn. 

If £(0) <0, let S,, be the unique solution of (S) + f£(0) =0, 
and S, the rheobase obtained as in the first case. If S < 8S, no excita- 
tion occurs; if 8. < S < S, there will be a finite number of dis- 
charges; if S = 8., there will be infinitely many discharges; and if 
S > 8., then after at most a finite number of discharges a state of 
permanent excitation will result persisting as long as the constant 
stimulus S is applied. 

This last case of f(0) < 0 might seem to be physically impossible, 
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but if it be supposed that an intensity S,, of the stimulus is injurious, 
then the mechanism breaks down for such stimuli and one still has 
only a finite number of discharges from any “physiologically admis- 
sible” stimulus. 

It is to be noted that on this theory the frequency of the dis- 
charge under constant stimulation is a fixed characteristic of the nerve 
fiber, and is independent of the intensity S of stimulation. Only the 
number of discharges varies with S, and this is potentially infinite 
for fibers for which f(0) < 0, but limited in those for which f(0) > 
0. In order to obtain variation of the frequency with S by a two-fac- 
tor theory (or, in fact, by an n-factor theory) it is necessary to gen- 
eralize equations (1) to a form 


d 
== fi (1 — 1°» tp — 22° 8); 

(52) 
d 
2 = f.(a— 2°, — 2°, 8), 


where it is supposed that the functions f vanish with their three argu- 
ments %, — 71°, % — x.° and S. If the functions f are expanded in 
power series of the three arguments, the linear equations discussed 
above may be regarded as first-order approximations. 

Nevertheless, the linear theory in the periodic case yields, at 
least qualitatively, not only the possibility of repetitive discharges, but 
also the depressed state (relative refractoriness) and the ensuing su- 
pernormal phase. However for any quantitative checks it will be neces- 
sary to study the solutions (41) with A and B constant giving the re- 
covery course after withdrawal of the stimulus, and also the course of 
development of the two substances with intermittent stimulation. 


LITERATURE 


H. A. Blair, 1932. “On the intensity-time relations for stimulation by electric 
currents,” Jour. Gen. Physiol., 15, 709-729. 

A. V. Hill, 1936. “Excitation and accommodation in nerve”, Proc. Roy. Soc. Lon- 
don, 119 B, 305-355. 

A. V. Hill, 1936. “The strength-duration relations for electric excitations of me- 
dullated nerve”, Proc. Roy. Soc. London, 119 B. 440-453. 

F. Offner, 1937. “Excitation theories of Rashevsky and Hill,” Jour. Gen. Physiol., 
21, 89-91. 

N. Rashevsky, 1933. “Outline of a physico-mathematical theory of excitation and 
inhibition”, Protoplasma, 20, 42-56. 

N. Rashevsky, 1938. Mathematical Biophysics, The University of Chicago Press. 

G. Young, 1937. “Note on excitation theories”, Psychometrika, 2, 103-106. 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 1, NUMBER 3 

SEPTEMBER, 1939 


DIFFUSION IN COLLOIDAL MEDIA 
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When the molecules of a solute diffuse through a medium con- 
taining large colloidal particles, which absorb the diffusing mole- 
cules, the latter are transported in the diffusion flow not as free 
molecules, but as absorbtion compounds: solute + colloid. When the 
colloidal particle is much larger than the molecule of the solute, and 
has therefore a much smaller mobility, this results in a reduction of 
the apparent diffusion coefficient for the solute. The biological impli- 
cations of this are discussed. 


In the first approximation the diffusion coefficient of a substance 
may be determined by considerations of the viscous resistance offered 
to the moving molecule by the solvent. In this way the well known 
formula of Einstein for the diffusion coefficient is obtained. (Gye- 
mant 1925). According to it, the diffusion coefficient should be rough- 
ly inversely proportioned to the average linear size of the diffusing 
molecule. Inasmuch as the permeability of the cell membrane is large- 
ly determined by the diffusion coefficient (Reiner, 1937) we would ex- 
pect the permeabilities to depend approximately also on the size of 
the diffusing molecule. However, it is well known that different mole- 
cules, having approximately the same size, permeate the cell some- 
times at very different rates. 

The present paper discusses the possibility that the diffusion 
process may be influenced by the adsorbtion of the diffusing molecule 
by larger colloidal particles. If the diffusing molecule is transported 
in form of such an adsorption compound, then its diffusion coefficient 
will be considerably lowered, when the carrying colloidal particle is 
very large. Of two molecules of same size, but different chemical na- 
ture, one may be strongly adsorbed, the other not. Thus we may find 
large differences in the diffusion coefficient and in permeability. 

The effective diffusion coefficient in a colloidal solvent may be de- 
termined upon the assumption that the solute and the colloid establish 
an adsorption equilibrium, while both experience the usual resistance 
to motion from the sclvent medium (e.g., water). The adsorption is 
calculated by means of the law of mass action. This amounts to assum- 
ing that the adsorption is stoichiometric. In the simple case, a mole- 
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cule of the solute (S) is adsorbed upon an active group of the colloid 
molecule (C) to form one molecule of adsorption compound (A), thus: 


StC2A+U (1) 


U is the heat of formation of a molecule of A. Let the concentrations 
of free S,C, A be n;, %-, % molecules cm; then the total concen- 
trations of solute and colloid respectively are: 


N,=M%+%; Neo=M+ MN. (2) 


The total diffusion current, J , of S (no. molecules of solute cross- 
ing per sec. an area of 1 cm.” normal to the flow) is the sum of the 
currents J, and J, of bound (A) and free (S) particles. For, if the 
distribution of N, is non-uniform, the same will in general be true of 
n, and n,, even though N, is constant (as will be assumed in the fol- 
lowing). Thus there will be a gradient of m, even when there is no 
gradient of total colloid N,. Moreover, by Fick’s law, 

J, = —D); erad i; 57. = — DD; stad 
where D, and D, are the diffusion coefficients of the substances S and 
A. Hence: 


J=J,+J,=—DgradN,, 
where, by Eq. (2) and the constancy of N.: 
D=D,[1— (0/0 Ns) ] + Da(0 1/0 Ns) (3) 


is the effective diffusion coefficient of the solute* S. From Eq. (3) it 
appears that D is the average of D, and D,, with a weighting factor 
0 Na/0 N,. It will be shown below that 0 < 0 n,/d N;, < 1 for all N, 
Ande N es 


For Eq. (1), the law of mass action takes the form: 
Note = NsNc 5 (4) 


where N, is the reciprocal of the equilibrium constant, and has the 
dimensions of a concentration. It may be interpreted as the N ¢ value 


for which n, = n,, i.e., for which half the solute is adsorbed. It is 
given by the Arrhenius equation: 


log No = log Z, — U/kT , (5) 


where k is the Boltzmann constant, T is the absolute temperature, and 


*Tf the diffusion coefficients of stubstances C and A are not equal, N, will 
not remain constant as assumed above. As C is a colloid, it will have very 1 
molecules compared to those of S ; it is thus reasonable to suppose that, ri a first 
approximation, D, = D,. If this is not the case, the above analysis will require 
generalization. 
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Z, a constant. Expressing , and », in terms of N,, N., and ™ by 
means of Eq. (2), one finally obtains the equation for nm, as an im- 
plicit function of N, and N,: 


Nota = (N. — Ma) (Ne — Ma) (6) 


To solve this, it is convenient to introduce the dimensionless quan- 
tities: 


ee en y=n,/Ne ; Z2=N,/N-. (7) 
It then takes the form: 
zy = (x —y) 1— y), (8) 


and has the solution: 


p= ee dey (2-1 2)? = 42). (9) 


The other solution, corresponding to the positive sign of the radical, 
is excluded on physical grounds: it makes y everywhere larger than x 
for z > 0, while, from Eqs. (2) and (7), x must always exceed y , and 
z > 0. Eq. (9) represents a one-parameter family of curves in the 
x-y plane (see Fig. 1). Regardless of the value of z , the amount of S 
adsorbed per molecule of C , y is zero for « = 0, and rises monotonic- 
ally toy =1forz =o. 

The parameter z is proportional to N,, which in turn, varies ex- 
ponentially with the heat of formation of A. The interpretation of 
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3(D5+D.) 


xX _ 
FIGURE 2 


Fig. 1 thus follows readily. As U approaches o (adsorption extremely 
exothermic), z approaches zero, so that y tends to a line consisting of 
the straight line y = x for 0 < x < 1, and the straight line parallel 
to the x-axis, y = 1, for « = 1. That is to say: if the reaction is pro- 
nouncedly exothermic, and C is in excess of S,(x < 1), all of S will be 
adsorbed; in the presence of an excess of S,(% s 1), an amount of S 
exactly equal to that of the C present will be adsorbed. For finite 
U,(z > 0) the amount of adsorption is less than this limit, so that 
there is always some free solute for every % < oo. 
The first derivative, y, = 0n./d N,, is given by: 
Ye = 3{1 — (w@—1+4 2)/[(e@+1-+ 2)? — 42]'}. (10) 


This family of curves is represented in Fig. 2. The function falls mo- 
notonically from 1/(z-+ 1) at the origin to the asymptotic value zero 
for x = o. On the other hand, it does not always vary monotonically 
with z, as shown in Fig. 3. If one differentiates Eq. (10) with re- 
spect to z, and sets the result equal to zero, y, is found to possess a 
maximum with respect to z. The value of z at this maximum is con- 
nected with x by the simple relation: 


oe il le Oa (11) 


Thus, since z > 0, the maximum can exist only for > 1. Fora <1, 
the variation of y, with z is monotonic. Eq. (11) can also be written: 


Nv, =N,—N.. (11’) 


Thus N, is the excess of solute S over colloid C at the z-maximum. The 
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maximum is accordingly attained for smaller excess of solute, the 
more exothermic the reaction. 

The previous assertion that 0 < 0 n,/d N. < 1 can now be justi- 
fied from Eq. (10). For the criterion that the second term of yz shall 
be no larger in absolute value than unity is found to be simply z > 0, 
which we know to be always true. Consequently, the entire expression 
in braces cannot become either negative or larger than 2, from which 


10 De 
9 
\ 
8 
PY 
6 
y,5 x=0 3 (B+ Dg) ? 
me x= 0.5 
sat x=1.0 
2 Sie esas fe 
x= 2.0 
2. x=4.0 
: D. 
“ 1 2 3 4 5 2 
Z— 
FIGURE 3 


follows the statement to be proved. Thus the expression (3) for D is 
seen to be a true average of D, and D,, lying between them in magni- 
tude. The variation of D may thus be simply obtained from that of yz 
by altering the scale; for D is a linear function of y., and D = D, 
when y, = 0, D = D. when y, = 1. Accordingly, the scale on the 
right-hand side in Figs. 2 and 3 gives D. If Da < D, (which is vir- 
tually always the case), D increases whenever Yz decreases and con- 
versely. D thus in this case has a minimum wherever y- (vide supra) 
has a maximum. 

The variation of D with z suggests a possibility of accounting for 
the conclusions of Landahl (Landahl, 1939) on the diffusion coeffici- 
ent in unfertilized Arbacia eggs. From available data, he calculates 
that D for oxygen may be as much as 10° times as large as D for lactic 
acid. Since D is the mean of D, and D, , these would have to differ by 
the same factor of 10° if the change in D is to be accounted for en- 
tirely by the change in y, , since the latter lies always between 0 and 
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1. This is not too improbable an assumption, for the values of colloidal 
coefficients of self-diffusion may fall extremely low, depending upon 
the density and structure of the particular colloid. 

Let us suppose that we are comparing the two solutes at the same 
concentration, and that the density of the colloid has not varied, so 
that.x is the same for both cases. Suppose, as is reasonable, x <1. 
Then the further assumption necessary is that the value of U is quite 
large for such a solute as lactic acid, and very low for a non-ionic and 
non-polar solute of the type O,, CO,. Then, from Fig. 3 it follows 
that the value of D for oxygen would lie upon one of the upper curves, 
far to the right (large z, small U) ; shifting to the lactic acid, with x 
unchanged, amounts to moving up the curve to the left, arriving at 
successively smaller values of D, which approach D, as a limit. The 
closer together the oxygen and lactic acid points lie on a given curve, 
the greater must be the difference between D, and D, to produce the 
desired change of D. Furthermore, such a change is less likely to be 
found if « > 1, since then y, varies only between 0 and 4. 

Adequate data on the U values do not appear available, and so it 
is necessary to assume the relation mentioned in the preceding para- 
graph. From the nature of the two types of molecule, such a hypothe- 
sis appears not unreasonable, though it is impossible to predict it with 
certainty on theoretical grounds. The difference of U values need not 
be excessively large; for, as shown by Eq. (5), No varies exponentially 
with U , and so may change considerably for relatively small changes 
of U. 

It is also possible to consider keeping the solute fixed and chang- 
ing the colloid. For « < 1, D will decrease if the colloid is replaced 
by one of a more active type (larger U for a given solute). For x > 1, 
the same will. hold up to a certain degree of colloid reactivity, after 
which a more active colloidal medium will cause an increase in the 
value of D. 

In a similar fashion, the variation of D with x may be obtained 
from the right-hand scale of Fig. 2. From this, D is observed to in- 
crease with increasing solute concentration, the increment for a given 
added amount of solute being sharpest in the neighborhood of x = 1. 
A caution is necessary, however, if this graph is to be applied for high 
values of x. For the principal assumption underlying the mass action 
law in the simple form of Eq. (4) is that the solution shall be dilute 
with respect to all constituents. Since x is a relative concentration 
by Eq. (7), it follows that the formulas are applicable for very high 
values of x only if N, is so low that, for the # under consideration, N 
still satisfies this restriction. al 
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